/* ------------------------------------------------------------------
		   Support code for problem set #4

			     SOURCE CODE
------------------------------------------------------------------- */

#include "support.hh"
#include <stdlib.h>
#include <sys/types.h>
#include <fstream.h>

matrix read_data(string filename)
{
  ifstream data(filename);

  int patterns;

  data.scan("# patterns = %d\n", &patterns);

  matrix M(patterns, 3);
  data >> M;

  return M;
}

int seeded = 0;

double random_number()
{
  if (!seeded) {
    srandom(getpid());
    seeded = 1;
  }

  return ((double) (random() % 100000)) / 50000.0 - 1.0;
}